1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.01
## LDevsum 1.00 1.00
## dh0 1.05 1.13
## dh1 1.06 1.17
## dl0 1.15 1.42
## dl1 1.38 1.91
## dl2 1.12 1.33
## dl3 1.12 1.33
##
## Multivariate psrf
##
## 1.3
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1447.5305 | 22552.6600 |
| DIC3 | 1304.7740 | 22449.1670 |
| PWAIC | 117.9956 | 136.0804 |
| WAIC | 1395.5854 | 22455.2355 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.01
## LDevsum 1.00 1.01
## dh0 1.94 3.64
## dh1 1.88 3.50
## dl0 1.27 1.71
## dl1 1.13 1.34
## dl2 1.08 1.22
## dl3 1.04 1.11
##
## Multivariate psrf
##
## 1.66
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1374.94287 | 27387.8319 |
| DIC3 | 1251.73037 | 27458.0214 |
| PWAIC | 90.74929 | 118.7216 |
| WAIC | 1320.99190 | 27462.1164 |